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Abstract 


A  rapid-compression  machine  is  a  laboratory 
apparatus  to  study  aspects  of  the  compression 
stroke,  combustion  event,  and  expansion  stroke  of 
an  Otto  cycle.  As  a  simple  model  of  such  a  machine, 
unsteady  one-dimensional  nonisobaric  laminar  flame 
propagation  through  a  combustible  premixture,  en¬ 
closed  in  a  variable  volume,  is  examined  in  the  as¬ 
ymptotic  limit  of  Arrhenius  activation  temperature 
large  relative  to  the  conventional  adiabatic  flame 
temperature.  In  this  limit,  a  thin  propagating 
flame  separates  nondlffusive  expanses  of  burned  and 
unburned  gas.  The  pressure  through  the  enclosure 
is  spatially  homogeneous  for  smooth  flame  propaga¬ 
tion.  However,  expansion  of  the  hot  burned  gas  re¬ 
sults  in  congressional  preheating  of  the  remaining 
unburned  gas,  and  in  fact  the  spatially  homogeneous 
gas  may  undergo  autoconversion  prior  to  arrival  of 
the  propagating  flame.  If  such  an  explosion  is  too 
rapid  for  acoustic  adjustment,  large  spatial  differ¬ 
ences  In  pressure  arise  and  the  resulting  nonlinear 
waves  produce  audible  knock.  Here  attention  is  con¬ 
centrated  on  what  fraction  (if  any)  of  the  total 
charge  may  undergo  autoconversion  for  a  given  oper¬ 
ating  condition,  and  what  enhanced  heat  transfer 
from  the  end  gas  would  preclude  autoconversion-- 
though  too  great  heat  transfer  from  the  end  gas 
could  result  in  flame  quenching  (unburned  residual 
fuel). 

K 

1.  Introduction 

During  nonisobaric  flame  propagation  through 
an  initially  homogeneous  premixture  of  combustible 
gases  in  an  (in  general,  variable-volume)  enclosure, 
end-gas  knock  may  occur.  That  is,  the  residual  un- 
bumed  gas  is  preheated  prior  to  flame  arrival  by 
compression,  owing  to  expansion  of  the  already 
burned  gas  and  possibly  also  owing  to  movement  of 
the  walls  of  the  enclosure;  if  this  preheating  re¬ 
sults  in  a  high-enough  unburned-gas  temperature, 
virtually  the  entire  residual  gas  may  undergo  auto- 
ignition  and  rapid  conversion  to  product  in  a  spa¬ 
tially  homogeneous  explosion.  This  phenomenon  is 
known  as  end-gas  knock;  because  the  temporal  scale 
of  the  abrupt  local  autoconversion  may  be  shorter 
than  the  acoustic-wave  speed,  the  pressure  field 
(well  approximated  to  be  a  function  of  time  only 
during  normal  flame  propagation)  becomes  spatially 
nonuniform,  and  the  resulting  nonlinear  waves  inter¬ 
act  with  the  walls  to  produce  loud  noise.  In  view 
of  the  current  capacity  to  treat  the  exhaust  of 
Otto-cycle  engines,  the  potential  damage  to  engine 
components  from  end-gas  knock  Is  a  major  obstacle 
to  operation  at  high  thermal  efficiency,  i.e. ,  at 
high  compression  ratio1**. 
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Here  further  insight  into  end-gas  knock  is 
sought  from  a  mechanical-engineering  point  of  view, 
as  distinct  from  a  chemical-kinetic  point  of  view. 
That  is,  a  Shvab-Zeldovich-type*  formulation  of  the 
one-dimensional  unsteady  laminar  interplay  of  spe¬ 
cies  and  heat  diffusion,  exothermic  chemical  reac¬ 
tion,  compressional  heating,  and  heat  loss  is  solved 
for  an  impervious  noncatalytic  variable-volume  en¬ 
closure.  Such  a  model  may  suffice  for  a  so-called 
rapid-compression  machine  for  cases  in  which  igni¬ 
tion  is  arranged  such  that  a  planar  flame  propagates 
parallel  to  one  moveable  (piston-like)  end  wall. 
Often,  such  a  machine  serves  as  a  conveniently  in¬ 
strumented  laboratory  apparatus  for  detailed  exami¬ 
nation  of  compression  and  reaction  (and  perhaps  sub¬ 
sequent  expansion  as  well).  In  contrast,  in  an 
Otto-cycle  engine  multiple  repetitions  of  the  en¬ 
tire  cycle  (including  induction  and  exhaust)  are 
entailed,  the  flame  propagation  is  highly  unsteady 
("turbulent"),  the  flame  is  often  predominantly  per¬ 
pendicular  to  the  piston  crown,  and  the  range  of 
parametric  variation  may  in  practice  be  constrained 
relatively  narrowly  by  cost  and  available  equip¬ 
ment1  »  11  . 

Accordingly,  the  following  questions  are  ad¬ 
dressed:  (1)  under  what  conditions  is  the  onset  of 
knock  to  be  expected  in  a  rapid-compression  machine; 
(2)  what  fraction  of  the  Initial  premixture  under¬ 
goes  autoconversion;  and  (3)  what  heat  transfer  from 
the  end  gas  would  suffice  for  smooth  flame  propaga¬ 
tion  across  the  entire  combustible  mixture. 

2.  Formulation 

For  a  low- Mach- number  flow  in  which  fuel  F  and 
oxygen  0  pass  exothermically  to  product  P  (in  the 
presence  of  nitrogen  N)  via 

Vp  F  +  Vg  0  +  vN  N  vp  P  +  vN  N  ,  (2.1) 

where  vi  is  the  stoichiometric  coefficient  of  spe¬ 
cies  i,  the  dimensional  conservation  equations  for 
mass,  species,  and  energy,  and  state  equation  are 
taken,  in  the  domain  0  <  x*<  L*(t*),  t*  >  0,  to  be5> 


l£l  ♦  =o-^=o*  - 

at7  3x*  u  ■  ax*  p  *  at* 

f9Y<  3Y<i  i  r  3Y-i 

°*[^»-  +  u*  ax7]  '  '5xT  3x*j  * 


-  p*u*  ;  (2.2) 
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In  particular,  a  direct  one-step  irreversible  bimolec 
ular  second-order  Arrhenius-type  global  reaction  is 
adopted,  and  no  attempt  is  made  to  examine  detailed 
chemical-kinetic  rates  and  mechanisms. 
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til* 

•  0  =  p  *  l  *  ; 


(2.12) 


p*  =  p*  R*  T*  .  (2.5) 

Here,  t*  and  x*  are  the  time  and  Cartesian  spatial 
coordinates,  respectively;  u*  is  the  gas  speed  and 
0*  is  the  streamfunction;  p*  and  T*  are  the  density 
and  the  temperature;  R*  is  the  gas  constant  for  the 
mixture,  taken  to  be  composed  of  species  of  compara¬ 
ble  molecular  weight;  p*  is  the  pressure,  taken  to 
be  a  function  of  time  only  from  consideration  of  the 
equation  of  conservation  of  momentum;  Y  and  m$*  are 
the  mass  fraction  and  molecular  weight,  respective¬ 
ly,  of  species  i;  m*  =  vr  flip*  +  vQ  mQ*  +  vr  m^;  D* 
and  X*  are  the  mass-transfer  and  thermal-conductivity 
coefficients,  respectively,  and  cp*  is  the  (constant 
universal'  heat  capacity  at  constant  pressure;  Q*  is 
the  heat  of  combustion  per  mass  of  premixture;  and 
w*  Is  the  reaction  rate.  The  Lewis-Semenov  number 
Is  taken  as  a  constant  of  order  unity,  i.e., 

(X*/c  *) 

Le  =  "(•p*b«)  =  const-  *  °(1)*  (2.6) 

It  is  also  taken  that 

p*D*  ^  p*n =  fnc(t*)  .  (2.7) 

Further,  consistent  with  the  nature  of  the  formula¬ 
tion. 


*  *°  *vF+vO  VF  v0  (  *  *  *  ,1  .  . 

w*  =  B ‘  T  p  F  °Yf  FYq  °exp|-Ta/(T  -Tuo)j,  (2.8) 

ir  ★ 

where  B1  is  the  (constant)  frequency  factor;  Ta  is 
the  (constant)  Arrhenius  activation  temperature;  Tj0 
Is  the  uniform  initial  premixture  temperature.  Intro¬ 
duced  to  modify  the  Arrhenius  factor  so  as  to  pre¬ 
clude  the  "cold-boundary  difficulty";  and  a  charac¬ 
terizes  the  preexponential  thermal  dependence  of  the 
reaction  rate.  Here,  a  =  {(vF+vQ)-l}  Is  adopted. 

It  is  convenient  to  adopt  the  von  Mises  trans¬ 
formation  (x*,t*)  *  (0*,t*),  where  the  streamfunc¬ 
tion  i F*  is  defined  in  (2.2).  Under  this  transforma¬ 
tion,  the  species  and  energy  conservation  equations, 
(2.3)  and  (2.4),  become 

3Y .  3Y .  *  «  Vp  vn 

3F  -  <p*  °*>  5*4  =  *  B<  YF  Y0  x 

expj-T*/(T*-TJo)),1=F,0,  (2.9) 
f£-  Le(p*V)  $£*  (^J  B‘*(pVR*)a  YfVfY0V°  x 

«p(-V/(T--T^)}  •  ;  (MO) 


Y  (vp ,  t  )  =  =  Up-  -2-U - 1  , 

YFu  ♦  Y0b 

T* ( tp* ,  t* )  -  T  * 

T(»p,t)  =  — T t-  .  y* -  ; 

bo  uo 

Pit)  •  , 

„(,)  -  edftgl ,  .  .“•Ifigl  , 

UO 


x(0,t)  = 


(2.13) 


(2.14) 


(2.15) 


Here,  the  subscript  u  denotes  unburned-premixture  con¬ 
ditions,  the  subscript  b  denotes  burned-state  condi¬ 
tions  at  the  flame,  and  the  subscript  o  denotes  ini¬ 
tial  conditions.  Specifically,  ujt,  is  the  initial 
adiabatic-flame  speed,  and  Tb*0  is  the  initial  adia¬ 
batic  flame  temperature.  Also  <f>  =  ( Y pu/yqu )  is  the 
equivalence  ratio,  with  0  <  l  for  cases  of  interest, 
so  Ypb=  0.  Hence,  with  0  -  (0/(1  -  0))  of  order  unity, 

Y0b  =  (V Ou  -  YFu>  =  an*d  Tb=  {TS  +  (Q*/CP*  >  YFu> 

(including  Tb„  =  Tu0  +  (Q*/cp  )Ypu  ).u  Note  also  that 

Po  =  Puo  P  ”uO’ 

Thus,  the  nondimensional  species  and  energy  equa¬ 
tions  are,  for  a,  n  =  1,  and  =  1, 

ay  ^  „  a2Y 
T r  ' EpFp 

=  -  -i-A0PYU  +  «)  exp{-  6(1  -  T)/T)  ,  (2.16) 

JT  -  e  Le  p  2fX 
TT  E  Le  p  FF 

--  j-  A0pY(l  ♦  *Y)  exp{-B(l-T)/T}  +  ^-Y^-^p-^, 

whprp  (2.17) 
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6  -  Ta*/(Tbo  -  TuV  >>l'  K  -  (Tbo  *  Vo^Vo*0^- 
Ao=  (B'*Duo/uuV)(po/R*,0t*"V°YFua  exp("B)  ” 1*(2-19> 


Y  ■  cp*/cv*  =  cp*/(cp* 


R*)  >  1  . 


(2.20) 


The  nondimensional  state,  continuity,  and  mapping 
equations  are 


„ .  p  au  _  a  1  .  a  f(i+KT) 

p " T+KT  ’  5^  '  [p  J  '  at  [  p  J  • 


ax-  1  -U^T) 

p  p 


(2.21a ,b,c) 


From  steady  (isobaric)  fuel-lean  laminar  flame 
propagation,  for  6»  1.  it  is  known  that,  for  \>F  =  1, 


The  mass  conservation  and  mapping  equations  are 


3u*  3  f  1  1  3x*  ,  1 

30*  Ft*  (p*-)  ’  1)0*  ~pr  ’ 

For  nondimensional ization, 


(2.11) 


a0  ■  rfr  \}  *  for  vf  " 1  •  {z-z2) 

where  the  higher-order  terms  are  known  but  are  not 
utilized  here. 

In  this  formulation,  the  fixed  wall  is  located 
at  x  *  0;  the  moving  wall  is  located  at  x  =  L(t),  with 
1(0)  =  1.  Under  the  von  Mises  transformation,  with 
the  fixed  wall  located  at  0  *  0,  the  moving  wall  is 


located  at  <i  -  t ( t) ,  with  y(0)  =  y0  =  const.  >  0, 
i.e.. 


V  =  I  p  dx.  =  p 


tj/  =  p  dx,  =  p 


(1+KT)  * 


iwry  • 


(2.23a) 


(2.23b) 


In  general,  the  initial  conditions  for  this  flow 
geometry  are 

Y  -*•  1,  T-*0  (p+1)  as  t-*-0  for  \ p  >  i|<fo-*-0  ,  (2.24a) 

Y-»0,  T -*•  1  (p-»l)  as  t-*-0  as  •  (2.24b) 

Note  that,  subject  to  (2.24),  it  follows  from  (2.23) 
that  Y0  =  1 .  In  general,  if  both  walls  are  imper¬ 
meable  noncatalytic  adiabatic  ones,  the  boundary 
conditions  are 

«  -yj— *0,  x,u-*0  as  ip+0  (t  >  0),  (2.25a) 

4 t-  ,  -gj-^O,  x-»  L(t) ,  u  +  l'(t)  as  tp -►  1  ( t  >  0) . 

**  ^  (2.25b)# 

The  jump  conditions  at  the  flames  are 


H  =  <VYU)--1.  [Tj-(Tb-Tu)- 

at  i)i  =  i))f(t). 


For  t>t-j,  the  unburned  region  consists  of  just  the 
end-gas  region.  To  distinguish  among  these  regions, 
the  following  notation  is  introduced: 

unburned  bulk  gas 

Y(i).t)  =  Yu(t|i,t) ,  T(i|/,t)  =  Tu(<)i,t}, 

Y  =  k  for  ( t)  <  i|/<  <|'i;  (3.3a) 

unburned  end  gas 

Y(«>.t)  =  Xu ( ip , t ) ,  T(i|>, t)  =(H>u(d-.t). 

Y  =  p  for  i)-j  <  i)  <  1 ;  (3.3b) 

burned  bulk  gas 

V(*.t)-Yb(*,t),  T(*,t)»Tb(«.,t), 

Y=  o  for  i)<i  >  i)if(t)  >  i)<  >  0;  (3.3c) 

burned  end  gas 

Y(i)st)  =  Xb(D,t),  T(Ui,t)  =  ®b(-i-,t), 

Y  =  T  for  <pf( t)  >  C'  >  .  (3.3d) 

For  the  unburned  bulk  gas,  (3.1)  and  (3.2)  become 
3Y  3T  ,  (1+K  T  )  . 

at  “  °*  3t  ‘  ic  Kp  It  °-  (3.4a,b) 

The  solutions  of  (3.4),  subject  to  the  initial  condi¬ 
tions  for  this  region, 

P(0)  -  1,  Y  (<p,0)  --  1,  Tu(*,0)  =  0,  and  (3.5) 

«Pf (0) 

with  tpf ( 0 )  =  it'fQ-*- 0,  ^  const,  (to  be  specified),  are 
.  .  .  .  l  17  .  o  (>c-1)/k  1 


(3.4a,b) 


3.  External  Hydrodynamics  of  Flame  Propagation 

Exterior  to  the  thin-flame  structure,  there  is 
an  unburned  region,  defined  by  <|/f(t)  <  <|»<  1,  and  a 
burned  region,  defined  by  i|if{t)  >  i)  >  0.  In  both 
these  unburned  and  burned  regions,  the  diffusion  and 
reaction-rate  contributions  are  taken  to  be  negli¬ 
gible,  so  that  the  species  and  energy  equations, 
(2.16)  and  (2.17),  reduce  to 

(3.u, 

Here,  the  unburned  region,  i|<f(t)  <  i) <  1,  is  considered 
to  consist  of  a  bulk-gas  region,  if^U)  <  and  an 

end-gas  region,  )|<t)i<l,  where  (given)  denotes 
the  Interface,  or  contact  surface,  between  the  bulk 
gas  and  the  end  gas.  This  division  of  the  unburned 
region  holds  for  0<t<tf.  When  t=tj,  all  the  un¬ 
burned  bulk  gas  has  been  burned,  i.e.,  <)f( t^ )  =  ijij . 

M 

On  physical  grounds,  since  pressure  is  spatially  in¬ 
variant,  results  depend  on  wall  sepa ration  as  a 
function  of  time;  as  long  as  the  temporal  history  of 
wall  separation  is  the  same,  results  are  Invariant 
to  which  (or  both)  end  wall  moves.  Also,  from 
(2.19)  and  (2.22),  altering  some  parameters  not  only 
alters  dimensionless  groups,  but  also  alters  the  de¬ 
rived  quantity  uUq,  so  that  the  temporal  scale  and 
the  streamfunction  scale  used  for  nondimens lonali- 
zation  are  altered.  Thus,  for  example,  altering  the 
value  of  the  activation  temperature  T-*  changes  the 
value  of  8  and  uu*0,  so  not  only  may  the  value  of  T 
at  any  value  of  ip  and  t  (in  general)  be  altered, 
but  the  associated  dimensional  quantities  ip*  and  t* 
are  altered.  Finally,  in  a  straightforward  general¬ 
ization  of  (2.24),  sparking  need  not  occur  neces¬ 
sarily  at  time  t  *  0;  e.g.,  congressional  heating 
owing  to  wall  motion  may  occur  for  a  finite  time 
Interval  before  sparking  initiates  flame  propaga¬ 
tion. 


1  17  1  (<-l)/K 

(2.26)  Yu(ip,t)  =  1,  Tu(ip,t)  =Tu(t)  = — ^ —  |p(t)j  -1  . 


In  a  similar  manner,  for  the  unburned  end  /j  6a  b\ 
gas,  the  solutions  of  ’ 

»V*  S®u  (u-1)  dp  . 

at  Tt - V  — Kp -  it  0  *  (3.7a,b) 

subject  to  the  initial  conditions 

p(0)  «  1,  Xu(ip,0)  *  1,  (5)u(ip,0)  =  0,  and  i)i  <  ip  <  1,(3. 8) 

are 

_  _  1  fc  s (w-l)/u 

M*»t)-i,  ©u(*.t)-  ®u(t)  =  -f- |p(t)J  -i  . 

(3.9a ,b) 

For  the  burned-bulk-gas  region,  (3.1)  and  (3.2)  become 
3Yh  3T.  .  (1+K  Th)  . 

ir=°- 17-^— kp — if10-  <3-10a*b> 

The  solutions  are  of  the  forms 

i  T  I  y (CT-l)/a 

Yb(».t)-Hb(*).  Tb(*,t)-4-Fb(*)|p(t)f  -1  . 


At  the  flame,  (3.11a, b) 

V*f(t),t)-Hb(*f(t))  :  Ybf ( t)  =  Hbf(t)  ,  (3.12a) 

i  r  i  s(o-D/o 

Tb(i|'f(t),t)  * Fb(i|/f(t)|p(t)j  -1 

Tbf(t)---r[Fbf(t){p(t)}{°  1)/°-1  •  (3.12b) 

For  t  <  t^  and/or  ipf ( t)  <i^,  the  jump  conditions  at  the 
flame  (between  the  unburned  bulk  gas  and  the  burned 
bulk  gas)  yield 


Ybf(t)  =  0  -»  Yb(ip,t)  =  0 


(3.13a) 


SR 
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3. 


f  T  wMff' '  iffliZl  *|  1121# 
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i  •  1+K  T  '  f  ■  l+K<Hl(4*.t)  ■ 

1  "pity  J- d*  ♦  L-:kiT--J  d» : 


pl>fU> 

L(t)=l - Wr  Fb(t)  diji  + 

■  p(t)  ■  ° 


(1-fc) 

V~p/+'ri^  •  ai5) 

Upon  multiplication  by  {p(t)}^°,  the  time  derivative 
of  (3.15),  with  rearrangement,  produces 


1  p-(t)  ,  (,-<■)  jvjH 

a  pit)  <  r  TI/ 

L(t)|p(t)| 


(U-O)  +  L 1  ( t )  .  KV(t) 


Ut){p(t)} 


Thus,  the  flame  speed,  the  difference  between  the 
speed  of  the  flame  front  and  that  of  the  unburned 
gas  at  the  flame  front,  is 

(  il  ^f'(t) 

S  «  S(t)  =  Xf'(t)  -  lim  <{u  (i^.t)  }■  =- - . 

(3.21) 

In  Appendix  A,  it  is  determined  that,  for  a,n  =  1  and 
vQ,Vp=  1,  the  flame  speed  is 

S(t)=|l +KTu(t)}|l  +  Tu(t)J2  exp  BTu(t)/2{l+Tu(t)} 

(3.22; 

It  may  be  shown  that  {3xu(i^,t)/3t}  =  u  (i^.t)  is 
continuous  at  i(i *<(*<,  i.e.,  is  continuous  across  the 
contact  surface  between  the  unburned-bulk-gas  region 
and  the  unburned-end-gas  region. 

Thus,  with  the  motion  of  the  moving  wall,  i.e., 
L(t)  and  L'(t)  and  the  Interface  location  tj^  speci¬ 
fied,  the  pertinent  initial-value  problem  becomes 
(3.6b),  (3.9b),  (3.16),  (3.21),  and  (3.22),  subject 
to 


U  LU 


L(t)p(t) 


=  0  .  (3.16) 


p(0)  =  1,  Tu(0)  =  0  , 


(3.23) 


Integration  of  the  mapping  equation  over  the 
unbumed  domain  <|i<^<1,  with  i|«>ij/p(t),  yields 

'L(t)  [*1  {l  +  KT  (V/j.t)} 

dx“  J - pltl - ,dV 


From  (3.16)  and  (3.23)  for  L'  =  0,  tc  =  o,  =  1,  one 
recovers  the  following  classical  results,  if  p=p, 
when  ipf  -  1 : 
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p(t)  -  1  =  oKi^p(t)  ■»  pt  -  1  -  oK  ■>  =  d'f(t) . 

This  implies  that  the  peak  value  of  Ty  is 
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This  problem  terminates  when  t  =  tj ,  such  that 

P(V  =  Pi’  Tu(ti>=Tui*  ®u(ti)  =  (H)ui^f(t1)^i. 

(3.24) 

where  p^,  Tu^ ,  (h)u^  are  found  in  the  course  of  so¬ 
lution.  here,  the  motion  of  the  moving  wall  is  ap¬ 
proximated  as 


L(t)  = 


./Ml  -  l\ 
2(CR) 


(1  -  COS  u  t)  , 


L'(t)  ° sin  in  t)  ,  ( 3.25a  ,b) 

where  ( CR) ( >1 )  is  the  compression  ratio  (L^x/t-min) 

=  (1/L,nin),  and  u  (=  2nfl  =  2tt{«*/ ( uj50/L0*  ) ) ,  with  ft* 
the  number  of  revolutions  per  unit  time)  is  the 
(nondimensional )  frequency. 

3.2  t )  >  vpj  (Flame  in  the  End  Gas) 


For  t>tj  and/or  tpf(t)>lPi>  integration  of  the 
mapping  equation,  (2.21a),  over  the  domain  0<<|i*l 
yields 

1  f"1 

L{t)=" - ~I7a  Fb(l(')d'F  +  7 - nyr 

J0  |p(t)} 


|p(t)| 


1 


■»l»f(t) 

Gfa(ij/)dil/  + 
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■  i-’Mt) ' 

- iW-  • 

{p(t)} 

Thus,  for  this  case. 


_i_£44 

o  Pit) 


(3.26) 
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♦f(t) 
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L  ( t )  ■ 

!p(t|} 

TTy 

«f(t) 


Gb(ip)d4,  -i 


L(t)|p(t)j  J 

where  it  is  approximated  that 


*V(t) 


(T -  a)  m  .  L'(t)  "Vf 

—a - 1 - TT7t  TltT  ‘  IXt)p7t)  '  °  * 


(3.27) 


fMt}  {«Pf(t)  ipi}  r  i 

Gb(q;)d^  =  ^ - 11  |Gb(4'f(t))+G>J(4/i)J  ,(3.28a) 

li 

(u-l)/u 


■  p(t)  ■ 


(u-l)/u 


J1 


Gb(¥<f(t))  = -1--  j^t-1)/T  ’  Gb^i*  *  1T-1)/T  • 

P(t)  p1 

v  >  (3.28b) 

Integration  of  the  mapping  equation  over  the  (end- 
gas)  unburned  domain  <  1 ,  with  \j»>i)if(t),  yields 

xu(4',t)  *  L(t)  -  ,  xu(^f(t),t)  =  xf(t)  * 

M 

{i  - 

no  -  r~w  •  (3-29a*b) 


In  turn, 

pp  (ij-.t)  =  uu(i|/,t)  =  L '  ( t)  + 


3t 


1  (1-0) 
p  r  il/u 


dx 


dt 


f  (t)  =  xf'(t)=L'(t)*^ 


{p(t!}_ 


(  l1/y 

jp(t)J 


«f'(t) 


,1/P 


{p(t)}1 
Thus,  the  flame  speed  is 
S  =  S(t) 


(3.30) 


xf'(t)  -  1  im 


{uu(U-,t)} 


t  +  Vflt) 

where,  for  a,n  =  1  and  vQ,Vp  =  1  , 
S(t)  =  {l  +  K©u(t)}{l+®u(t)}2 

exp  6®u(t)/2{l  +  <H)u(t)} 


V  (t) 


’J  {P(t,} 


\l/v  * 
3.31) 


(3.32) 


Thus,  for  this  domain,  the  pertinent  initial- 
value  problem  is  given  by  (3.9b),  (3.27),  (3.28) 
(3.31),  and  (3.32),  where 

P(ti)  =  pi,®u(ti)=  ®ui,  .  (3.33) 

The  problem  nominally  terminates  when  t=ti,  such 
that 

P(tj)  «  pr  (jDuU,)  =  ®ul.  =  1  .  (3.34) 

The  above  generalization  of  earlier,  variable- 
volume-enclosure  work7  entails  a  flame  speed  which 
is  derived  from  the  adopted  reaction  rate  (cf.  Sena- 
chin  and  Babkin8)  and  which  reduces  to  the  laminar 
isobaric  value  as  t-*0  (cf.  Shivashinsky*) . 

4.  Discussion 


First,  in  Appendix  B  it  is  shown  that  for  cir¬ 
cumstances  of  Interest,  heat  transfer  (to  slippery 
noncatalytic  side  walls)  may  be  simulated  conveni¬ 
ently  by  assignment  of  values  to  the  polytropic  con¬ 
stants  k,  p ,  a,  and  T;  i.e.,  the  polytropic  constants 
effectively  are  equivalent  to  introduction  of  heat- 
transfer  coefficients.  The  smaller  the  value  as¬ 
signed  to  a  polytropic  constant,  the  greater  the  heat 
transfer  In  the  associated  domain.  While  adequate 
experimental  identification  of  polytropic  constants 
(or  heat-transfer  coefficients)  seems  yet  to  be  re¬ 
ported  in  the  literature,  the  consequences  of  assign¬ 
ing  relative  values  permits  comparison  of  heat- 
transfer  requirements  for  varying  parametric  condi¬ 
tions. 

Second,  in  singular-perturbation  terms,  end-gas 
knock  is  the  physical  manifestation  of  the  possible 
nonuniform  validity  In  time  of  the  asymptotic  ap¬ 
proximation  (based  on  B»l)  that  the  reaction-rate 
contribution  is  negligible  in  the  unburned  gas;  how¬ 
ever,  the  asymptotic  approximation  (based  on  c«  1) 
that  the  diffusion  contribution  is  negligible  remains 
uniformly  valid.  Thus,  if  one  retains  chemical  re¬ 
action  terms  in  (2.16)  that  were  discarded  in  writ¬ 
ing  (3.1,2)  because  the  terms  were  taken  as  expo¬ 
nentially  small,  then  one  has 


(4.1) 


if’  -  ~c~  Ao  pY(1  +  *v>  expj -6(  1-T)/t|  , 

TF=~T  AopY  (1  +  ,”Y)  exp|-B(l-T)/T}  + 

(7-1)  (1+KT)  dp 
7  Kp  dt  * 


(4.2) 


where  (Y(4,t),  T(4,t),  y)  -  (Yu(t),  Tu(t),  k)  in  the 
unburned^bulk-gas  region  and  (r(<l\t),  T(4,t),  7 )  ♦ 
(Xu(t),(0u(t),  u)  in  the  unburned-end-gas  region. 

It  is  recalled  that  adjustment  of  the  polytropic  con¬ 
stant  7  simulates  the  effect  of  heat  loss.  The 


above  problem,  subject  to  the  initial  conditions 
Y{ 0)  =  1,  T ( 0)  »  0,  is  a  variant  of  the  classical  spa¬ 
tially  homogeneous,  nonstationary  ignition  problem, 
in  which  the  criterion  for  an  explosion  selfsus- 
tained  against  heat  loss  is  sought.  Here  congres¬ 
sional  heating,  as  reflected  in  the  pressure  p(t) 
raises  the  temperature  of  the  premixture;  clearly,  as 
T  rises  in  time,  the  reaction-rate  term  may  become 
comparable  to  the  other  terms  in  the  energy  equation. 
If  one  substitutes  the  pressure,  p(t),  from  solution 
of  the  problem  posed  in  Section  3,  then  one  can  ob¬ 
tain  an  indication  from  (4.1)  and  (4.2)  whether  sig¬ 
nificant  autoconversion  occurs  in  the  end  gas  prior 
to  completion  of  flame  propagation.  If  modest 
autoconversion  causes  pressure  rise  In  time  dur¬ 
ing  flame  transit  beyond  that  accounted  for  in 
the  theory  of  Section  3’,  that  development  is  not 
easily  generalized  to  include  it.1*  Of  more  concern 
is  very  significant  autoconversion  of  unburned  gas 
on  a  time  scale  comparable  to,  or  less  than,  the  time 


scale  for  completion  of  flame  propagation.  Espe¬ 
cially  if  the  time  scale  is  less  than  that  for  com¬ 
pletion  of  flame  propagation  (and  less  than  that 
for  acoustic  adjustment),  there  occurs  effectively 
a  constant- volume  explosion  of  the  unburned  gas. 

Such  a  relatively  Instantaneous  reaction  of  unburned 
gas  results  in  a  pressure  in  the  end-gas  region 
higher  than  the  pressure  in  the  rest  of  the  gas. 

The  consequence  Is  nonlinear  pressure  waves,  and  an 
audible  sound  from  sympathetic  vibration  of  engine 
components.  However,  there  is  no  Incentive  to  pur¬ 
sue  solution  Into  this  spatial ly-nonuniform-pressure 
domain  because  the  above  formulation  is  Inadequate. 


Standard  asymptotic  solution  for  B»1  of  (4.1) 
and  (4.2),  or  even  of  (4.2)  with  reactant  consump¬ 
tion  Ignored  (i.e.,  Y  =  1 ) ,  suffices  to  establish 
that,  for  plausible  parameter  assignments  for 
hydrocarbon-air  premixtures  of  automotive  interest, 
autoconversion  becomes  effectively  instantaneous 
relative  to  flame- propagation  speed  over  (say)  a 
half  centimeter  for  unburned  gas  temperature  on  the 
order  of  1025  K;"  practical  experimental  data  for 
nonlsobaric  combustion  also  give  this  value.1  Mathe¬ 
matical  treatment  beyond  that  in  Appendix  C  seems 
of  limited  engineering  utility. 


Thus,  the  following  computational  procedure  is 
evolved  to  answer  the  questions  posed  in  the  last 
paragraph  of  Section  1.  First,  w’th  ic.o  taken  as 
given  and  <|u  set  equal  to  unity  (so  there  is  no  end 
gas),  the  first  initial-value  problem  posed  in  Sec¬ 
tion  3  Is  carried  out  for  the  parametric  assign¬ 
ments  of  Interest.  If  It  hold  uniformly  in  time 
that  the  temperature  of  the  unburned  gas  TJ<  1025  K 
(for  specificity),  then  knock  Is  taken  not  to  occur. 
If,  however,  Tj  reaches  1025  K  for  ijif(t)  <  1,  then 
that  value  of  Is  taken  as  <.<;  I.e.,  Is  the 
fraction  of  the  charge  Identified  as  bulk  gas  and 
(l-^j)  as  end  gas.  Then,  with  T  taken  as  given,  by 
trial  afld  error,  the  largest  value  of  u  Is  sought  for 
which  ©  -»  1025  K  as  4if*l.  That  is,  the  minimal 


amount  of  heat  transfer  from  the  end  gas  is  sought 
such  that  significant  end-gas  autoconversion  barely 
is  precluded.  In  this  manner,  the  fraction  (if  any) 
of  the  charge  that  would  undergo  rapid  autoconversion 
is  identified,  and  the  relative  heat-transfer  provi¬ 
sion  to  preclude  such  autoconversion  is  determined. 

In  the  present  formulation,  with  adiabatic  end 
walls,  heat  transfer  is  effected  through  (noncataly- 
tic  slippery)  side  walls.  Incorporation  of  finite 
heat  transfer  through  end  walls  would  necessitate 
introduction  of  thermally  diffusive  layers  in  the 
vicinity  of  ip  =  0,1;  while  such  boundary  layers  would 
thicken  in  time,  they  would  remain  thin  relative  to 
the  flow-domain  length,  for  the  time  span  and  the 
parametric  assignments  of  interest  (e«l).  Such 
end-wall  quench  layers  have  been  examined  in  other 
contexts  ;  their  introduction  in  the  present  context 
would  seem  to  incur  nonessential  complication,  such 
as  end  gas  in  a  nonuniform  thermodynamic  state.  The 
only  precaution  (already  noted)  is  that  very  exces¬ 
sive  heat  transfer  from  the  end-gas  could  so  retard 
flame  propagation  that  expansional  cooling  (from 
piston  motion)  would  leave  residual  unburned  pre¬ 
mixture  on  time  scales  of  practical  interest. 


Computational  Examples 


A  limited  number  of  examples  now  are  presented 
of  the  procedure  discussed  in  the  penultimate  para¬ 
graph  of  Section  4.  For  brevity  of  reference,  the 
following  assignment  of  (mainly)  dimensional  param¬ 
eters  is  termed  nominal: 

T*q  =  300  K,  T*q  =  1800  K,  T*  =  15,000  K,  <t>  =  0.8, 

Duo  =  0,1  cn,2/s-  %0  =  40  cm/s*  Lo  =  10  cm*  Le  =  !• 

ft*  =20  rev/s.  (5.1) 

It  is  consistent  with  these  assignments  to  adopt  the 
following  "nominal"  values  for  dimensionless  param¬ 
eters  (the  parameter  c  does  not  enter  the  calcula¬ 
tions— its  value  is  given  only  for  completeness): 

8=10,  A  =  50,  K  =  5,  4  =  4,  e  =  2«10'\  (0*31.4, 


y  =  1.4,  k=  1.4,  T=  1.4,  a  =  1.4,  CR=  6,  0.  =  135°. 

19  (5.2) 

The  last  assignment,  alluded  to  earlier  in  the  foot¬ 

note  to  (2.25b),  implicitly  gives  the  time  of  spark¬ 
ing,  more  precisely,  the  time  at  which  flame  propa¬ 
gation  commences.  That  is,  it  is  defined  that 

6  =  ut,  (5.3) 


so  at  the  beginning  of  compression  t  =  0,  it  holds 
that  6  =  0,  which  corresponds  to  "bottom  dead  center" 
[L  =  1,  the  maximum  value  of  L,  from  (3.25a)].  At  6  = 
0ja,  assigned,  i.e.,  at  t,g  =  0jg/(o,  the  variables 
\|)At)and  x*(t),  held  at  zero  for  Oct*  tig,  take  on 
values  consistent  with  (3.21)  and  (3.18).  While 
sparking  and  compression  could  commence  simultaneous¬ 
ly  so  tig  =  0fg  =  O,  in  general  sparking  occurs  after 
a  finite  Interval  of  compression,  so  t^g>0.  In  the 
nominal  case,  6^g  =  135°,  or  flame  propagation  com¬ 
mences  after  about  (6/7)  of  the  length  of  the  (com- 
pressional)  stroke.  All  parameters  are  held  at 
their  nominal  values  through  the  computations  unless 
it  is  explicitly  stated  otherwise. 

Onset  of  knock  is  taken  to  occur  at  Tu  =  0.4, 
or,  from  (5.1),  at  a  temperature  of  900  K--a  value 
slightly  lower  than  that  cited  elsewhere  in  this 
paper. 


6 


The  nominal  case  is  found  (by  design)  to  be 
just  barely  knock-free.  Computation  indicates  that 
combustion  is  completed  at  ^  0.088;  since  the  basis 
of  nondimensionalization  is  (Lg/u*0)  = 0,25  s,  com¬ 
bustion  lasts  22  ms,  or  about  three-to-four  times 
longer  than  in  a  typical  Otto-cycle- type  automotive 
engine.  At  the  completion  of  burning,  TF 0.399, 
p  -  46.6,  e  =  157°  ,  L  =  0. 197,  whereas  at  the  start  of 
burning  t  =  0.075,  Tu  =  0 . 129 ,  p  =  5.69,  L  =  0.289.  As 
with  all  cases  to  be  examined,  combustion  is  com¬ 
pleted  during  the  compression  stroke,  i.e.,  the 
crank  angle  9  <180°;  of  course,  in  an  automotive 
context  the  combustion  interval  is  likely  to  be 
(very  roughly)  symmetric  with  respect  to  the  time  of 
piston  position  at  top  center. 


If  the  polytropic  constant  of  the  burned  gas  o 
is  reduced  in  value  from  1.4  to  1.1,  all  other 
parameters  being  held  fixed,  then  at  the  completion 
of  burning  t  =  0.089,  =  0.378,  p  =  40.945,  0  =  159.8°, 

L  =  0.192,  whereas  at  the  start  of  burning  the  cor¬ 
responding  values  are  (of  course)  those  given  in  the 


preceding  paragraph.  Since  modifying  the  heat 
transfer  from  the  burned  gas  is  seen  to  have  but 
modest  effect  on  the  temperature  of  the  unburned 
gas,  this  parametric  variation  is  not  pursued  here 
further. 


Table  1  presents  results  for  four  parametric 
variations,  all  selected  to  incur  knock.  The  heat- 
transfer  requirement  from  the  end  gas  that  pre¬ 
cludes  knock  is  identified  in  the  form  of  the  maxi¬ 
mum  value  of  u  (denoted  pcrit;)  that  still  permits 
(H}j  <  0.4  at  all  times  during  the  combustion  event. 
Further  results  for  the  case  CR =  12  are  presented 
in  Figs.  1-3,  and  further  results  for  the  case 
K=3,  in  Fig.  4-6;  the  other  two  variations  exam¬ 
ined  are  0j_q  *  165°  and  w=  51.6.  These  cases  entail 
between  13?  and  57%  of  the  initial  combustible  mass 
being  involved  in  autoconversion  (in  the  absence  of 
augmented  end-gas  heat  transfer),  and  in  this  sense 
represent  a  spectrum  of  cases.  On  Figs.  2,  3,  5, 
and  6  it  may  be  recalled  that  x;(t),  the  interface 
between  the  unburned  end  gas  and  the  unburned  bulk 
gas,  is  not  defined  for  Xf(t)>x^(t),  and  accord¬ 
ingly  is  neither  computed  nor  plotted.  It  may  be 
noted  that  the  larger  the  fraction  of  end  gas  (i.e., 
the  smaller  the  value  of  <(/j),  the  larger  the  heat- 
transfer  requirement  to  preclude  knock  (i.e.,  the 
smaller  the  value  of  Ucrjt)-  However,  it  is  note¬ 
worthy  that  reduction  of  K  (from  its  nominal  value 
of  5)  to  3,  while  8  is  held  constant,  implies  T*  is 
altered  from  300  K  to  450  K;  since  Tu  =  0.4  is  re¬ 
tained  as  a  knock  criterion,  T* =  990  K  becomes  the 
temperature  for  explosive  conversion  of  reactant  to 
product  gas. 

As  a  cautionary  note  concerning  other  paramet¬ 
ric  variations,  it  may  be  noted,  for  example,  that 
if  Ypu  is  altered,  then  $  and  t£0  are  altered,  from 
the  discussion  below  (2.15);  hence,  8,  A,  and  K  are 
altered  from  (2.19)  and  (2.22),  but  also  uj„  (uti¬ 
lized  in  the  nondimensionalization  of  time  t*  and 
the  piston-speed  parameter  O*)  is  altered  fr\.m 
(2.19). 

6.  Concluding  Remarks 

The  present  development  has  concentrated  on  en¬ 
hancing  heat  transfer  as  a  countermeasure  to  end- 
gas  knock.  Another  (perhaps  complementary)  strategy 
is  to  seek  fuel  additives  that  inhibit  end-gas  ajto- 
conversion  without  impeding  flame  propagation. 

While  trial -and-error  testing  yielded  the  knock- 
inhibiting  fuel  additives  of  the  past  (such  as  tet¬ 
raethyl  lead,  now  constrained  from  use  because  it 


poisons  exhaust-gas-treating  catalysts),  it  remains 
to  be  seen  whether  retention11  of  detailed  chemical- 
kinetic  rates  and  mechanisms  for  fuels  of  practical 
interest  can  assist  identification  of  a  substitute, 
environmentally  acceptable  knock-inhibiting  addi¬ 
tive.  ** 


Appendix  A.  Flame-Zone  Structure 


As  in  the  unbounded  (isobaric)  case,  the 
structure  of  the  flame  zone  for  the  bounded  (non- 
isobaric,  enclosed)  case  for  e«l,  with  8»1, 
consists  of  two  regions:  a  convective-diffusive 
region,  and  a  diffusive-reactive  region.  The 
results  presented  are  for  a,n  =  1  and  vp,v0  =  1. 


For  the  convective-diffusive  region,  the  ap¬ 
propriate  transformation  of  variables  is  ( *p , t )  •+ 
U,t),  where 


5U,t;...)  =  M(t)E(t;8) 


4f(t)-  v)/cj  • 


( A.  1  ) 


The  domain  of  this  region  is  -“  <  4  <  0,  t>0,  such 
that  the  unburned  exterior  region  is  approached  as 
Further,  the  flame  speed  is  taken  to  be 


S(t; . . . )  -  '4f’(  t;  . . )/ (p(t))1/K  « 

V(t)E(t,8)[l  +  0(8'ir|  •  (A. 2) 

Here, 

E(t;8)  =  exp|e{Tb(t)  -  l}/2Tb(t)l  = 


expl8Tu(t)/2(l  +  Tu(t))  ; 


(A. 3) 


M(t)  =  (Tb(t))2  =  U  +  Tu(t)}2,  (A. 4a) 

V(t)  =  {p(t))(lc'1)/,c{Tb(t))J  ={l+KTu(t)Hl+Tu(t)}J. 

(A. 4b) 

Note  that  C "*■  ( ->^)/e  and  S-*l  as  t->-0,  such  that  the 
unbounded  case  is  the  initial  condition  of  the 
bounded  case.  It  is  through  the  matching  of  the 
solutions  of  this  convective-diffusive  region  with 
those  of  the  diffusive-reactive  region  that  it  is 


The  evolution  of  a  planar  flame  propagating 
through  a  not-too-far-from-stoichiometric  premix¬ 
ture  enclosed  in  a  long-enough  constant-volume  duct 
to  a  so-called  tulip12’13  shape  may  be  evidence  of 
a  "mini-knock"  event.14  It  is  observed  that  such  a 
flame  ultimately  propagates  less  rapidly  in  a  duct 
at  midheight  (relative  to  near-ceiling  and  near¬ 
floor  flame  position).  Although  the  temperature  in 
the  congressional ly-heated  unburned  gas  may  not 
exceed  about  600  K  in  such  tests,  perhaps  cool- 
flame  (and  possibly  other  subsequent,  modestly  exo¬ 
thermic,  low-temperature)  reactions  engender  enough 
end-gas  counterflow  to  the  deflagration-wave  speed 
that  flame  progression  is  retarded  at  duct  midheight; 
presumably  friction  inhibits  the  counterflow  near 
ceiling  and  floor,  so  the  flame  locally  propagates 
relatively  unimpeded.  Incidentally,  for  long-enough 
tubes,  or  for  tuDes  of  small -enough  cross  section, 
spatial  gradients  of  pressure  plausibly  arise  such 
that  the  tulip-flame-shape  effect  could  be  observed 
even  if  the  extremity  of  toe  tube  near  the  burned 
(or  even  unburned)  premixture  were  open  to  the  atmos¬ 
phere.12  Since  the  end-gas  temperatures  are  typi¬ 
cally  over  1000  K  for  knock,  the  significance  of 
cool-flame  (and  possibly  other  subsequent  low-tem¬ 
perature)  reactions  in  rapid- compress  ion  machines 
under  knocking  conditions  is  unclear. 


demonstrated  that  (A. 2)  is  the  correct  asymptotic 
representation  of  the  flame  speed. 

For  this  region,  the  dependent  variables  are  of 

the  forms 

Y(<|>,t;...)  =  xU,t;...)£l+0(e'‘)]  ,  (A. 5a) 


T(*.t;...)  =Tb(t)  -  6'1{Tb(t)}2G(n,t;...)[l  +  0(fr1)j 


k,t;...)  =  FU,t;...)[l  +  0(fT1)j 


(A. 5b) 


Thus,  for  these  scalings,  the  specific  time-deriva¬ 
tive  terms  and  the  reaction-rate  terms  can  be  ne¬ 
glected,  such  that,  to  leading  order  of  approxima¬ 
tion,  (2.16)  and  (2.17)  become 

0-  ir=  0  *  (A-6) 

Le@--Sr=°-  ^ 

The  boundary  conditions  for  these  equations  are 
X 1 ,  F-»Tu  as  C -*•-<»  ;  (A. 8a) 

X-+0,  F-*Tk=  (1  +  T„)  as  C-0.  (A. 8b) 


X 1 ,  F-»Tu  as  C -*•-«>  ;  (A. 8a) 

X-*0,  F- Tb  =  (1  +  Tu)  as  5-0.  (A. 8b) 

The  solutions  of  these  convective-diffusive-region 
boundary-value  problems  are 


X(£,t)  =  ^  -  exp^fj, 

F(£,t)  =  jju(t)  +  exp(5/Le)j  • 

Tfa(t)  -  £  -  expU/LeTJ  . 
It  is  seen  that 


X  -*■  1,  F-*-Tu  (exponentially)  as  ?-►-« 


and  that 

Xn,-  5(  1  +...)-*■  0+  , 

Fn,Tb+(e/Le)(l  +  ...)-Tb  as  6 +■ 0- 


(A. 10) 

(A. 11) 

(A. 12a) 
(A. 12b) 


For  the  diffusive-reactive  region,  the  trans¬ 
formation  of  variables  is  (i/»,t)  * (n,t),  where 

n(*.t,...)-E(t;3)|{*f(t)-*}/B-lc]  ,  (A. 13) 

where  E(t;f!)  is  given  in  (A. 3).  The  domain  of  this 
region  Is  -«» <  n  <  °°  ,  t>0,  such  that  the  burned 
exterior  region  is  approached  as  n-+".  The  flame 
speed  Is  still  that  given  in  (A. 2)  et  seg.  Again, 
note  that,  as  t  +  0,  n-+  ( xe  ancTS^T. 

For  this  region,  the  dependent  variables  are 
taken  to  be 

Y(i|i,t;. . .)  =  8~l(Tb(t)}2  Z(n,f,...)[l  +  0(e-‘jj  - 
6*‘{l  +  Tu(t)>J  Z(n,t;...)[l  +  0(8-‘)],  (A. 14a) 


=  a  +  Tu(t)  J-B-Hl  +  Tu(t)}2G(n,t;.  .  .  )x 


1  +  0(8-')  . 


(A. 14b) 


For  e«l  and  g»l,  under  these  scalings,  to 
leading  order  of  approximation,  for  this  reqion, 
the  convective  terms  can  be  neglected,  and  (2.16) 
and  (2.17)  become 


3^t  -  1  exP(-G). 

(A. 15) 

Le  w  -  2L?*  2  • 

(A. 16) 

The  boundary  conditions  are 

Z,G  ->0  as  n  ■+  ®  ; 

(A. 17a) 

Z, G  -*■  ®  as  n  . 

v A. 17b) 

By  subtraction  of  (A. 16)  from  (A. 15)  ant 
(A. 17),  Z  =  LeG.  Then  (A. 15)- (A. 17)  red 

of 

C 

W  =  G  exp(‘G)  ; 

(A. 18a) 

G  +  0  as  n  +  “,  G  +  »  as  r)  + 

(A. 18b) 

The  first  integral  of  this  equation  is 

(A. 19) 

Thus,  it  is  determined  that 

G  n.  (exp(-n/2  ^  Le)}{l  +. . . )  +  0  as  n  ■+  ®  ;  (A. 20a) 


G  -(n/Le)(l  +. . . )  as  n  -*■  -® 


(A. 20b) 


Upon  returning  to  the  original  dependent  vari¬ 
ables,  it  is  seen  that,  for  the  convective-diffusive 
region  and  diffusive-reactive-region  solutions 

Y  n,  |j-s)  +  ..rj  +...  as  t,  -  0-  , 

Y  n,  B"1  ju+Tu)2{(-n)  +  ...)  +...  as  n  +  .  (A. 21a) 

T  n,  1+TU)  -  (-C/le)+...j  +...  as  c  -  0-, 

T  n,  ( l+Tu)  -  6"1  f(l+Tu)2{(-n/Le)  +  . ..}]  +  ...  as 

L  J  (A. 21b) 

Since  n=  B£/(1+Tu)2,  the  solution  for  these  two 
regions,  which  involved  the  specification  of  the 
leading-order  approximation  of  the  flame  speed, 
match.  In  the  text,  this  flame  speed,  expressed  as 

S(t)  =  S(Tu(t);8,K)  = 

(l+KTu(t)Hl+Tu(t)}2  exp[8Tu(t)/2{l+Tu(t))],  (A.22) 
is  employed. 


Appendix  B_. _ Pol ytrop i c_Law  to  Simula te  JHeat 

Transfer  Puri",  one  Combustion  Event 
in  an  Enclosure 

If  the  use  of  a  super  asterisk  to  denote  a 
dimensional  quantity  is  dispensed  with  in  this 
appendix,  then,  in  conventional  notation,  from 
thermodynamics , 


dT 

"3t~ 


d ( 1/P )  _ 
dt 


(B .  1 ) 


where  q  denotes  the  rate  of  heat  addition  per  unit 
mass  in  the  enclosure.  Hence,  for  an  ideal  gas, 

d(kn  T)  ,  n  d(?n  p)  _ 
dt  '  Ky~1'  dt 


djtn  p)  _  d(£.n  p)  =  q 


dt 


dt 


cvT 


(B.2) 


be  neglected  in  the  end  gas,  although  diffusive 
transfer  remains  negligible.  This  statement  may  be 
more  readily  noted  by  rewriting  (2.16)  and  (2.17) 
under  (C. 1)-(C. 3)  as 

P  +  <TV )  exp{ -6£  (  ©t-  -  T)/T} ,  (C.4) 


il 

3t 


p  Y  ( 1  +  *Y)  exp{-Be(©E-  T)/T)  + 


Y-l  ( 1  +  KT )  dp 
y  Kp  3t  ‘ 


(C.5) 


For  0^T<  ©e,  the  reaction  term  is  negligible; 
this  interval  is  also  characterized  by  the  follow¬ 
ing  : 


0  <  t  <  t£  ,  l<p<n£  , 


(C .  6) 


For  tractabil ity,  let  q  =  -e  cv  T,  with  the  con¬ 
stant  e  >0;  the  product  cc  plays  the  role  of  a 
heat-transfer  coefficient. v  Then,  if  subscript  o 
denotes  the  initial  state, 


d 

It 


•  fin  [ — E — 

-  Ip^JJ 


£ 


-L  =  -2-  exp(-c  t) 
D  Do 


(B.  3) 


During  combustion,  as  a  rough  approximation,  if  the 
constant  p  >  0,  with  p  distinct  from  that  of  (3.3b), 


_E_ 

Po 


= exp(pt) 


(B.4) 


hence,  substitution  for  t  by  means  of  (B.4)  gives 


exp(-ct)  =  expH-  -~ 


*n[-f-ll 
1  Mo 


,  -  --E- 
J  Ip, 


-e/p 


Use  of  (B. 5)  in  (B.3)  gives 


JE_ 

P, 


l+(e/p) 


ILl)U=  Ijl) 


I/’T* 


■  V 


1 

y* 


l  + 


(B.5) 


(B.6) 


Since  (e/p)  >  0,  y>Y*-  That  is,  the  effect  of  heat 
transfer  from  a  flow  with  rising  pressure  may  be 
simulated  by  adoption  of  the  polytropic  constant 
reduced  from  the  adiabatic  value.  Furthermore, 
via  (B.6),  the  polytropic  constant  may  be  related 
to  the  heat- transfer  coefficient  (and  a  characteri¬ 
zation  of  the  rate  of  pressure  rise). 

Appendix  C.  End-Gas-Autoconversion  Criterion 

In  (2. 16) -(2.22) ,  with  e  «  1  and  6  »  1 ,  it  is 
anticipated  that  a  physically  interesting,  auto¬ 
conversion-related  case  entails  the  following 
quantitative  relation  between  the  values  of  the 
(more  generally,  independent)  parameters  c,B: 

exp(-BT£),  with  T£  =  0(1).  (C.l) 

It  is  convenient  to  introduce  the  notation 

®eE(l  +  Te)-‘<l^TeMl-©e)/©e  ; 

6C  =  (1  +  Te)8>6  .  ( C .  2 ) 

Hence,  from  (C.l)  and  (C.2), 

e  =  6e  exp{-8£  (1  -  ®£))  (C.3) 

The  subscript  e  denotes  conditions  at  which  it  is 
anticipated  that  chemical  conversion  may  no  longer 


where  the  definitions  of  the  constants  te  and  1^  is 
evident.  Within  this  interval,  solution  is  given 
by  ( 3 . 6a , b )  in  the  unburned  bulk  gas,  and  by 
(3.9a,b)  in  the  unburned  end  qas.  For  the  typical 
values  e  =  2xl0"‘' ,  8  =  10,  (c.l)  gives  Te  =  1.31. 

(HV  =  0.432,  Be  =  23.1;  for  e  =  2xl0-4  but  6=5  (a  some¬ 
what  small  value),  TE  =  2.35,  H  =  0.299,  BE  =  16.8. 
From  (2.13,  if  Tj0  =  300  K,  Tbo  =  2100  K,  the  condi¬ 
tion  T<H^  implies  T*  =  10^8  K  for  6=10,  T  =840  K 
for  8=5.  The  value  for  T*  used  in  the  text  is  on 
the  order  of(H^  =0.4,  or  T*  =  1025  K. 

Equations  (C.4)  and  (C.5)  may  be  brought  into 
the  classical  canonical  form  for  study  of  transition 
to  explosion  by  introduction  of  the  following  expan¬ 
sion  about  the  state  e: 


t  =  B_ 


(i+*)(ne/®E) 

_ 


-  T_ 


,  y  =  i-b£-,©£x 


t  =  H  E  (i  +  6£-1®)  ,  p  =  nc  ( i  +  8£ - 1  n). 

Substitution  of  (1.7)  and  (C.8)  into  (C.4)  and 
(C.5)  yields  to  leading  order: 


(C.7) 

(C.8) 


4t'~  exp(@) 


3®  _ 

-yp-=exp(®)  + 


(Y-l)  (UK®c)  dn 
Y  K®e  'ar 


where 
X  =  ® 


hzl) 

y 


(l  +  K®  ) 


n 


(C.9a,b) 


(C. 10) 


is  an  integral.  To  lowest  order  at  the  outset  of 
the  transition,  reactant  depletion  need  not  be  ac¬ 
counted  for  in  the  heat  balance.  Also,  it  is  re¬ 
called  that  adjustment  of  the  value  of  y  is  taken 
as  a  convenient  alternative  to  the  introduction  of 
a  heat-transfer  coefficient  for  incorporating 
thermal-loss  effects  during  a  process  in  which  the 
pressure  rises  in  time.5§ 

T5 - 

Compressional  heating  exceeds  heat  loss  to  walls 
in  the  end  gas,  for  cases  of  automotive  interest. 
Hence,  the  end  gas  is  supercritical  in  the  sense 
that  pressure  and  temperature  are  assured  to  rise 
in  time.  The  concern  is  whether  or  not  the  in¬ 
crease  of  temperature  and  pressure  in  time  (from 


9 


Incidentally,  from  (C.l),  (C.2),  and  (C.8),  for  a  '"Johnson,  J.  H.,  Myers,  P.  S.,  and  Uyehara, 
long-enough  (say,  adiabatic,  constant-volume)  duct,  0.  A.,  "End-gas  temperature,  pressures,  reaction 
£•♦0,  Te  increases  and(ft  decreases;  in  such  cases,  rates  and  knock,"  SAE  Technical  Paper  650505, 

chemical  conversion  could  supplement  significantly  Society  of  Automotive  Engineers,  Warrendale,  PA, 

compression  as  a  mechanism  for  unburned-gas  preheating.  1965. 
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the  compressional-preheating/heat-loss  balance) 
is  such  that  transition  to  explosion  is  deferred 
so  that  completion  of  flame  propagation  occurs 
instead.  Thus,  while  representing  heat  loss  by 
adjustment  of  the  polytropic  constant  precludes 
study  of  the  near-critical  state  in  which  chemi¬ 
cal  exothermicity  and  heat  loss  balance,  the 
present  formulation  is  not  aimed  at  such  already 
well -studied  questions. 


Fig.  1  For  nominal  parametric  assignments,  defined 
by  (5.1)-(5.3),  except  that  the  compression  ratio 
CR= 12 ,  the  dimensionless  unburned-bul k-gas  tempera¬ 
ture  Tu  and  the  dimensionless  unburned-end-gas  tem¬ 
perature  are  presented  as  functions  of  dimen¬ 
sionless  time  t,  where  compression  commences  at 
t  =  0.  When  the  flame  enters  the  end  gas,  Tu  is  no 
longer  defined.  Here,  the  scale  for  nondimens ion- 
al ization  of  time  is  22  ms;  T  =  0  corresponds  to 
300  K,  T=0.2,  to  600  K,  T  =  0.4  to  900  K.  Augmented 
heat  transfer  from  an  appropriate  portion  of  the 
mixture,  the  last-to-burn  charge  (end  gas),  pre¬ 
cludes  any  unburned  gas  achieving  a  rapid-autocon- 
vers ion-precipitating  temperature,  taken  as  900  K. 
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Fig.  2  For  the  case  of  Fig.  1,  the  dimensionless 
dependent  variables  pressure  p,  combustion-chamber 
length  L,  and  interface  position  Xj  between  un¬ 
burned  bulk  gas  and  unburned  end  gas  are  presented 
as  f'  ictions  of  dimensionless  time  t.  When  the 
flame  enters  the  end  gas,  x^  is  no  longer  defined. 
Pressure  is  nondimens ionali zed  against  its  initial 
value,  here  1  atm,  and  lengths  are  nondimensional- 
ized  against  the  initial  chamber  length,  here  10  cm. 


Fig.  4  For  the  nominal  parametric  assignments  de¬ 
fined  by  (5.1)-(5.3),  except  that  the  ratio  of 
chemical  exothermicity  to  initial  enthalpy  K  =  3, 
the  dimensionless  temperature  Tu  and  the  dimension¬ 
less  unburned  end-gas  temperature^!^  are  presented 
as  functions  of  time  t.  The  scale  for  nondimension 
alization  of  time  is  22  ms;  T=0  corresponds  to 
300  K,  T  =  0.2  to  620  K,  T  =  0.4  to  090  K  (which  is 
taken  as  the  rapid-autoconversion-temperature  to  be 
avoided  by  augmented  heat  transfer  from  the  end 
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Fig.  3  For  the  case  of  Fig.  1,  the  dimensionless 
dependent  variables  chamber- length  L,  interface  po¬ 
sition  xj  between  unburned  bulk  gas  and  unburned  end 
gas,  the  flame  position  Xf,  and  the  fraction  <j*t  of 
total  initial  premixture  burned,  are  presented  as 
functions  of  time  t.  The  length  xf,  'ike  xit  is 
nondimens ionali zed  against  initial  chamber  length, 
here  10  cm. 


Fig.  5  For  the  case  of  Fig.  4,  a  presentation 
exactly  analogous  to  that  of  Fig.  2. 


